Chiral loop corrections to strong decays of positive and negative parity charmed 

mesons 
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We consider chiral loop corrections to the effective couplings which describe the strong interactions 
of two heavy mesons of even or odd parities and a light pseudoscalar meson, within a framework 
which combines heavy quark effective theory and chiral perturbation theory. The strong couplings 
are extracted from the experimental results for the decay widths, after including chiral corrections. 
We find that the introduction of positive parity heavy meson fields into the framework gives relevant 
contributions to the chiral loop corrections of the effective couphngs. The counterterms' contribu- 
tions are studied by making a randomized fit distribution of the relevant couplings at /i ~ 1 GeV. 
Then we discuss the chiral extrapolation based on calculated loop contributions. The inclusion of 
positive parity mesons in the chiral loops is found to be important for the lattice QCD studies of 
charm meson properties, where the knowledge of chiral limit behavior is required. 

PACS numbers: 13.25.Ft, 12.39. Fo, 12.39.Hg 



I. INTRODUCTION 

In recent years, discoveries of open charm hadrons have 
stimulated many studies. The first evidence of broad 
states in the charmed spectrum was provided by the 
CLEO collaboration [l[ which observed a broad £)^(2460) 
state with the features of an axial meson. In 2003, 
Belle and Focus Q reported the observation of fur- 
ther broad resonances D*+ and Dl°, ca. 400-500 MeV 
higher above the usual D states and with opposite par- 
ity. In the same year BaBar Q announced a narrow 
meson Dsj{2i\l)'^ . This was confirmed by Focus Q 
and CLEO [g| which also noticed another narrow state, 
Dsj(2463)+. Both states were also confirmed by Belle 
The measured properties of the Dq and D'^ states support 
their interpretation as belonging to the (0+, 1+) multiplet 
of cu and cd mesons. Conversely, the D^j states have 
been proposed as members of the (0+, 1+) spin-parity 
doublet of cs mesons [1, Q . 

Studies of the basic properties of these states have been 
triggered particularly by the fact that the £)sj(2317)+ 
and Z?sj(2463)+ states' masses are below threshold for 
the decay into ground state charmed mesons and kaons, 
as suggested by quark model studies El and lattice 
calculations [13, The strong and electromagnetic 

transitions of these new states have alr eady been stud- 
ied within a variety of approaches [1, [ij, [H, [H, [13, 
m, [H, [10, [m, [13, [H. In these investigations heavy 
hadron chiral perturbation theory (HHxPT) at leading 
order was used as well in attempts to explain the observed 
strong and electromagnetic decay rates [13, [IE, [lO]- In 
ref. [2J| the authors have considered the leading and sub- 
leading {l/mu) operators which describe the even-parity 
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charmed meson finite mass corrections within heavy me- 
son chiral perturbation theory. This calculation con- 
tained 11 unknown parameters in its predictions. After a 
number of attempts to determine these parameters, the 
authors of [13| found that a fit to the even-parity masses 
is possible, provided the couplings describing the strong 
interaction between even and odd parity heavy meson 
states and light pseudoscalars are constrained to lie be- 
tween and 1. 

In ref. [1^ the chiral loop corrections to the D* — > Dtt 
and D* —f Dj decays were calculated and a numerical 
extraction of the one-loop bare couplings was first per- 
formed. Since this calculation preceded the discovery of 
even-parity meson states, it did not involve loop contribu- 
tions containing the even-parity meson states. The ratios 
of the radiative and strong decay widths, and the isospin 
violating decay D* Dgir^ were used to extract the 
relevant couplings. However, since that time, the experi- 
mental situation has improved and therefore we consider 
in this paper the chiral loop contributions to the strong 
decays of the even and odd parity charmed meson states 
using HHxPT. In our calculation we consider the strong 
decay modes given in TABLE (T] The existing data on 
the decay widths enable us to constrain the leading or- 
der parameters: the D* Dtt coupling g, DqDtt coupling h, 
and the coupling g which enters in the interaction of even 
parity charmed mesons and the light pseudo-Goldstone 
bosons. Although the coupling g is not yet experimen- 
tally constrained, it moderatelly affects the decay ampli- 
tudes which we consider in this paper. 

Due to the divergences coming from the chiral loops 
one needs to include the appropriate counterterms. 
Therefore we construct a full operator basis of the rel- 
evant counterterms and include it into our effective the- 
ory Lagrangian. We first systematically calculate the chi- 
ral corrections to the wave function renormalization and 
the relevant vertex corrections. Then we constrain the 
three strong couplings, first by including the loop effects 
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Meson 




Mass [GeV] 


Width [GeV] 


Br. [%] (final states) 


D*+ [26] 


1" 


2.010 


(9.6 ±2.2) X lO^'^ 


67.7 ±0.5 iD°-K+), 30.7 ±0.5 (-D+vr") 


D*° [26] 


1" 


2.007 


< 0.0021 


61.9 ± 2.9 (D°7r°) 


[3] 


0+ 


2.403 ± 0.014 ± 0.035 


0.283 ± 0.024 ± 0.034 


(!)%+)" 




0+ 


2.350 ± 0.027'' 


0.262 ± 0.051*" 






1+ 


2.438 ± 0.030" 


0.329 ± 0.084^= 





"Observed channel 
''Average of Belle 
'^Average of Belle 



and Focus H values from 12711 
and CLEO [l| values from [27 



TABLE I: Experimentally measured properties of excited charmed mesons used in our calculations. 



while neglecting finite contributions from the countert- 
erms. The reason why one can expect their contribu- 
tions not to be significant is related to the appropriate 
choice of the renormalization scale and the fact that the 
initial and final meson states in processes we consider do 
not contain strange quarks. Still we study the impact of 
the counterterms by taking their values to be randomly 
distributed, and observing the results of the couphngs' 
fit. 

A full calculation of the strong decay couplings should 
contain, in addition to the corrections we determine, also 
the relevant I/itlh corrections as discussed in ref. j2^ . 
However, the number of unknown couplings is yet too 
high to be determined from the existing data. Also, re- 
cent lattice QCD studies [2§| of the strong couplings 
of heavy mesons have noticed that these corrections seem 
not to be very significant but pointed out the importance 
of controlling chiral loop corrections. 

Due to the increase of simulation time, when approach- 
ing the chiral limit, lattice studies use large values of 
the pion mass. In order to connect their results to the 
physical limit they employ the chiral extrapolation. The 
knowledge of chiral loop contributions is essential for the 
stability of this procedure. We compare two scenarios 
of chiral extrapolation within HHxPT, one including the 
contribution of the even-parity states, and the other with- 
out their contributions. 

The paper is organized into the following sections: Sec- 
tion II contains the description of the framework. Section 
III is devoted to calculations. In Section IV we explain 
the extraction of the bare couplings, while Section V con- 
tains the study of chiral extrapolations. We conclude 
with discussion in Section VI of the paper. 



II. FRAMEWORK 

We use the formalism of heavy meson chiral La- 
grangians [13, [3l|. The octet of light pseudoscalar 
mesons can be encoded into E = = cxp(2i7r*AY/) 
where the tt^A* matrix contains the pseudo-Goldstone 
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and / « 120 MeV at one loop |32|. The heavy- light 
mesons are customarily cataloged using the total angu- 
lar momentum of the light degrees of freedom in the 
heavy meson which is a good quantum number in 
the heavy quark limit due to heavy quark spin symme- 
try. The negative (j^) = l/2~ and positive [jf) = 1/2+ 
parity doublets can then be respectively represented by 
the fields H = 1/2(1 + y)[P;7'' - ^75], where P* and 
P annihilate the vector and pseudoscalar mesons, and 
S = 1/2(1 + y)[Pr^7''75 - Pq] for the axial-vector (P^*^) 
and scalar (Pq) mesons. 

The Lagrangian relevant for our study of chiral correc- 
tions to strong transition among heavy and light mesons 
is then at leading order in chiral and heavy quark expan- 
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>C;^ ±'^-1- -\- £'1+ + 'Cmix, 

-Tl- [Ha{lV ■ Vat - Sab^H)Hb] 

+gTT [HbHaAabJs] , 

Tr [Saiiv ■ Dab - Sab^s)Sb] 

±.gTr [SbSaAab^b] , 

hTl [HbSaAablb] + h.c. 



(2) 



ab 



5abd^ — is the covariant heavy meson deriva- 
tive. The light meson vector and axial currents are de- 
fined as = 1/2(^^9^^ + id^S,'^) and = i/2(e^9^e - 
^9^^^) respectively. A trace is taken over spin matrices 
and the repeated light quark flavor indices. Accordingly 
is of the order 0{p^) in the chiral power counting 
while the rest of this leading order Lagrangian is of the 
order 0{p^). Exceptions are the A// and A5 residual 
masses of the H and S fields respectively. In a theory 
with only H fields, one is free to set /S.h = since all 



3 



loop divergences are cancelled by 0{mq) counterterms 
at zero order in 1/mjy expansion. However, once S fields 
are added to the theory, another dimensionful quantity 
IS.SH = As — A// enters loop calculations and does not 
vanish in the chiral and heavy quark limit [24| . It is of the 
order 0{p'^) in the standard power counting and in order 
to keep a well behaved pcrturbativc expansion we need 
to modify this prescription formally and assign a power 
counting of l^sH ^ ■ At the end of the calculation we 
fix its value close to the phenomenological mass splitting 
between the even and odd parity heavy meson multiplets 
Asif ~ 400 MeV although smaller values have also been 
proposed when taking into account next to leading order 
terms in l/mn expansion [23 |. Conversely, when per- 
forming an extrapolation of light quark masses to the 
chiral limit, we need to treat these non-vanishing Ash 
contributions carefully. 

Following Refs. [2^, [s^, we absorb the infinite and 
scale dependent pieces from one loop amplitudes into the 
appropriate counterterms at order 0{mq) 
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Tr [HaHaAcl^{ml)c 
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-Tr [HaHbiv ■ PbcAaTs] 
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AiTr [SaSbim'^)ba] - A'^Tr [SaSaim^q)bb 



' A2 
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Tr [(55/75 )a6( 



■m\)ba 



g^3 
A2 
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+ -7-5- Tr [S'aS'aAc75(m|)cfc] + -r^Tr [S'cS'a(m^)ab Ac 75] 
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-Tr \SaShiv ■ Vbc^A c 



>75 



■a^ 

'AT 
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. -^Tr[(7?5/75)a6l 
Tr [{HSA-iz)aa{m\)bb\ 

TV [HaSaAcl5{ml),b] 



m^)ba 



and A-^ = 47r/. Ellipses denote terms contributing only 
to precesses with more then one pseudo-Goldstone bo- 
son as well as terms with (iv ■ V) ac ting on H or S, 
which do not contribute at this order [25| . The matrix 
ruq ~ diag(TO„, ?Tic;, 771s) induces masses of the pseudo- 
Goldstone mesons to^^ = 4Ao(777q + mb)/ f"^, where a, 5 
are the light quark flavor indices while m| = {(,mq£^ + 
S.^tTiq^^^). Thus, in our power counting (777^ ~ p"^) all 
the A terms in C^*^ arc of the order 0{p^) while the S 
and K terms are of the order 0{p^). Parameters X[ and 
X[ can be absorbed into the definition of heavy meson 
masses by a phase redefinition of H and S, while Ai 
and Ai split the masses of SU{3) flavor triplets of Ha 
and Sa, inducing residual mass terms in heavy meson 
propagators: A^ = 2Ai777a and A^ = 2Xima respec- 
tively [2^. As with A^fs), only differences between 
these 0{p^) residual mass terms enter our expressions. 
We denote them as Aba = Af, — A^, A^;^ = A^ — Aq 
and Af,a — A^ — A^ and fix them from phenomenolog- 
ical mass splittings between heavy mesons of different 
SU{3) flavor. Assuming exact isospin and heavy quark 
spin symmetry, the only nonvanishing splittings are then 
A3a « Ag^ « A35 ~ Asa « 100 MeV, where here 
a = 1,2. For the ki and Kg terms only the combina- 
tion KiQ = Ki + Kg will enter in an isospin conserving 
] manner here [25| (the ki — Kg combination contributes 
to isospin violating D* —> Dgir'^ decay, which we do not 
consider in this paper) . In the same manner wc only con- 
sider contributions of k[q = k[ + Kq and kig = ki + kg. 
At any fixed value of rUq , the finite parts of K3 , K3 and 
Tr [HaHbipbcV ■ Aa75] '^3 can be absorbed into the definitions of g, g and h 
respectively psl ]. However, one needs to keep in mind 
that these terms introduce a non-trivial mass dependence 
on the couplings when chiral extrapolation is considered. 
The 62 and ^3 enter in a fixed linear combination, in- 
troducing momentum dependence into the definition of 
g of the form g — {S2 + ^s)^ • k/A^. For decays with 
comparable outgoing pseudo-Goldstone energy, this con- 
tribution cannot be disentangled from that of g p5l | . On 
the other hand, these contributions have to be consid- 
ered when combining processes with different outgoing 
pseudo-Goldstone momenta. The same holds for contri- 
butions of 62 and 63 in respect to g, as well as S2 and 5'^ 
in respect to h. At order 0{mq) we are thus left with 
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bb\ 



— Tr [SaSbiTPbcV ■ A 
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explicit analytic contributions from 715, Kg, K5, Kig, k\ 
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^19; 



h-52 



and 82 



2 rp^ 

X 

h.c. - 



HaSbiv ■ VbcAcald] 



A2" 

X 

S' 

-^Tr 
Av 



Tr [HcSaim^q)abAbcl5] 



III. 



HaSbipbcV ■ Acalb] 



CALCULATION OF CHIRAL LOOP 
CORRECTIONS 



Here V^^^A^^ = 9^ Ac + [V'',^^]a6 is the covariant 
dcrivate acting on the pseudo-Goldstone meson fields 



Here we present the most important details of our chi- 
/rjNral loop calculations. We employ dimensional regulariza- 
tion in the renormalization scheme where the subtracted 
divergences 2/e — 7 -|- log47r -|- 1 are absorbed into the 
appropriate counterterms. 
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Wavefunction renormalization 



FIG. 2: Sunrise road topology diagram 



We first calculate the wave function renormalization respectively, which yield identically 



Z2H of the heavy H = P, P* and Pq, fields. This is 
done by calculating the heavy meson self-energy Il{v ■ p), 
where p is the residual heavy meson momentum, and 
using the prescription 



2H 



1 dn{v ■ p) 



2 dv ■ p 



on mass— shell 



(4) 



■^ab\a 
167r2/2 



3g"C — ,m. 



/ Aba — ^SH 



(6) 



where the mass-shell condition is different for the H and 
S fields due to their residual mass terms Ajj(5) as well 
as for the different SU{3) flavors due to the Aq terms. 
In general it evaluates to v ■ p — Ah(s) — Aa = 0. 

At the 0{p^) power counting order we get non-zero 
contributions to the heavy meson wavefunction renor- 
malization from the self energy (" sunrise" ) topology dia- 
grams in FIG. [1] with leading order couplings in the loop. 
In the case of the P mesons both vector P* and scalar 
Pq mesons can contribute in the loop yielding for the 
wavefunction renormalization coefficient 
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X' X' 

V m,. 



ba 



-,m, 



(5) 



As in Ref . [25| , a trace is assumed over the inner repeated 
index(es) (here b) throughout the text, while the loop 
functions Ci and their analytic properties are defined in 
the Appendix B. In our power counting scheme, their an- 
alytic Aba dependence is of the order 0{p'^) or higher and 
could be neglected at this order. However A^a also enter 
non-analytically and we have to check for sensitivity of 
our results to these parameters. On the other hand our 
power counting prescription for Ash ~ leads to well 
behaved 0{p^) and C(p^) analytic contributions to Ci. 
At leading order in heavy quark expansion, due to heavy 
quark spin symmetry, the wavefunction renormalization 
coefficient for the P* field is identical to that of P al- 
though it gets contributions from three different sunrise 
diagrams with states P, P* and Pj* in the loops [s^ . 

The positive parity Pq and Pj* obtain wavefunction 
renormalization contributions from self energy diagrams 
(FIG. [I]) with P*, P and Pq, Pj*, P* mesons in the loops 



B. Vertex corrections 

Next we calculate loop corrections for the PP*tt, 
PoPj^TT and PqPtt vertexes. At zeroth order in 1/toq 
expansion these are identical to the P*P*7r, P*Pj*7r and 
P*P*TT couplings respectively due to heavy quark spin 
symmetry. Again we define vertex renormalization fac- 
tors for on-shell initial and final heavy and light meson 
fields. Specifically, for the vertex correction amplitude 
T(v • Pi,v ■ pf,k^) with heavy meson residual momentum 
conservation condition pf = p^ + one can write the 
renormalization coefficient ZiHiHj-n schematically 



HiHfTT 



1 



T{v-pi,v-pf,kl) 



Ti,o.{v-p,,vpf,kl) 



on mass — shell 



(7) 



Here F/.o. is the tree level vertex amplitude, Pi(f) is 
the residual momentum of the initial (final) state heavy 
meson Hi(^f^, while k-^ is the pseudo-Goldstone mo- 
mentum. This implies that in the heavy quark limit, 
one must evaluate the above expression at fcj = mj. 



and consequently v ■ ky, 



Ai. where A 



«(/) 



is the residual mass of the ini- 



tial (final) heavy meson fields. While such prescription 
violates Lorentz mo mentum co nservation at the order 
Odfe^l/m//) w ©(y^A^- — ml/niH); where A:^ is the 

pion three-momentum, it ensures that in all expressions 
one only encounters the physical, reparametrization in- 
variant quantities Ash, Aab, Aab and Aab Satisfy- 
ing botli conditions would require the introduction of two 
separate velocity scales as in P ^ 13 meson transitions 
(see e.g. [1^), which is beyond the scope of this work. 

At the 0{p^) order, the relevant contributions to the 
vertexes between the heavy and light mesons come from 
the "sunrise road" loop topology diagrams in FIG. [2] with 
leading order couplings in the loop. Contributions from 
other one loop topology digrams either yield zero due 
to Lorentz invariance or are cancelled off in our further 
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calculations by the light meson wavcfunction rcnormal- 
ization Z^^i |3j| at the given order. For the case of the 
PP*-n vertex, only {P*,P), (P*,P*) and (Pq-A*) con- 
tribute pairwise in the loops. Adding the relevant 0{p^) 
counterterm contributions we thus obtain 



Zip, J 



\j \i \j 



^db 



-c 



A^b + A 



SH 



Kc{mq)cb,^^ , cab^. ^ Aba §2 + 83 



KKb 



K 9 



(8) 



The same expression is obtained for the P*P*tt vertex 
renormalization from pairs of {P,P*), (P*,P*), {P*,P) 
and (Pj* , Pj*) running in the loops. 

Similarly for the PoPi tt and PiP^tt vertexcs we get 
contributions from pairs of (Pi*,Po), (Pi*, Pi*), (P,P*) 
and (Po, Pr), (Pi*, Po), (Pi*, Pi*), (P*,P*) respectively 
running in the loops which yield identically 



, = 1 - 



KcKd^db 
2 f 2 



A„, A, 



g Ci , ,771 



-C 



As - Ash 




+ 



Kcimq)cb f^, ,s:ab~ X Aba ^2+^ 
T-— (,^.19 + K5j ; . 

KKb 9 



(9) 



Finally for the PPqtt and P*PiTt vertexes the pairs 
of (Po,P), (P*,Pi*) and (Pr,P*), (P*, Pi*), (P,Po) con- 
tribute in the loops respectively, yielding for the renor- 
malization identically 



= 1 - 



MicKd^db 



I o I Aca Adb 
X S ^99^1 , , 777^ 



Aca- Ash A^,+Ash 




Kcimq)cb 

KKb 



(I .x<^b/^ Aba - Ash S'^ + 5'^ 

(-19 +^5 .5) —■ 



(10) 



IV. EXTRACTION OF BARE COUPLINGS 

Using known experimental values for the decay widths 
of D+* , D^* , and , and the upper bound on the 
width of D°* one can extract the values for the bare 
couplings g, h and g from a fit to the data. The decay 
rates are namely given by 



r(p: ^ TT^n) 



\9f. 



6np 



■\K 



(11) 



and a similar expression (up to polarization averaging 
phase space factors) for r(Po — > ttP) and r(P* — > ttP*) 
with g coupling replaced by h and replaced by 

E^\k.^\. Here k-^ is the three-momentum vector of the 
outgoing pion and Et^ its energy. The renormalization 
condition for the couplings can be written as 



s. _ „ V^2P„v/^2P-\/^ _ 

9p'P,,n' - 9 / 9^ 



Zi 



9 

P*Pb^' 



(12) 



Pa P.* 



with similar expressions for the h and g couplings. 

Due to the large number of unknown counterterms en- 
tering our expressions (ks, K19, K5, K19, K5, k'^q, 62 + 5^, 
62+83 and 52+53) we cannot fix all of their values. There- 
fore we first perform a fit with a renormalization scale set 
to ^ ~ 1 GeV and we choose to neglect counterterm 
contributions altogether. Our choice of the renormaliza- 
tion scale in dimensional regularization is arbitrary and 
depends on the renormalization scheme. Therefore any 
quantitative estimate made with such a procedure can- 
not be considered meaningful without also thoroughly 
investigating counterterm, quark mass and scale depen- 
dencies. 

We constrain the range of the fitted bare couplings by 
using existing knowledge of their dressed values and as- 
suming the first order loop corrections to be moderate 
and thus also maintaining convergence of the perturba- 
tion series: 

• g - following quark model predictions for the pos- 
itive sign of this coupling [361 as well as previous 
determinations [3, [la, HI, l25l . [28l . [29j we constrain 
its bare value to the range g € [0, 1]. 

• h - this coupling only enters squared in our expres- 
sions for the deca y ra tes and was recently found to 
be quite large [T^.Tlsl. [l6l. [23j . We constrain its bare 
value to the region \h\ G [0, 1]. 

• g - non-relativistic quark models predict this cou- 
pling to be equally signed and smaller than g [STf . 
Similar results were also obtained using light-cone 
sum rules [l^, while the chiral partners HHxPT 
model predicts \g\ = \g\ d, H, HI]- A recent lattice 
QCD study [1^ found this coupling to be smaller 
and of opposite sign than g. We combine these dif- 
ferent predictions and constrain the bare g to the 
region g G [—1,1]. 

We perform a Monte-Carlo randomized least-squares fit 
for all the three couplings in the prescribed regions us- 
ing the experimental values for the decay rates from TA- 
BLE U to compute ^^'^ using values from PDG [2^ 
for the masses of final state heavy and light mesons. In 
the case of excited Dq and D[ mesons, we also assume 
saturation of the measured decay widths with the strong 
decay channels to ground state charmed mesons and pi- 
ous {D^ Dtt and D[ -> D*tt) [l^. We obtain the 
best-fitted values for the couplings g = 0.66, \h\ = 0.47 
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and g = —0.06 at = 3.9. The major contribution to 
the value of comes from the discrepancy between de- 
cay rates of D[ and Dq mesons. While the former favor a 
smaller value for \h\, the later, due to different kinemat- 
ics of the decay, prefer a larger \h\ with small changes 
also for the bare g and g couplings. Similarly, as noted 
in Rcf. (39j . such differences arc due to the uncertain- 
ties in the measured masses of the broad excited meson 
resonances. We expect these uncertainties to dominate 
our error estimates and have checked that they can shift 
our fitted values of the bare couplings by as much as 20% 
depending on which experimental mass values are consid- 
ered. On the other hand our results are fairly insensitive 
to the value of the residual mass splittings between heavy 
fields as 30% variations in Ash, Aaf,, Aab and Aab only 
shift our fitted values for the bare couplings a few per- 
cent. We have probed our results to the sensitivity to 
the renormalization scale /i and obtained a moderate de- 
pendence, namely a 20% variation of scale around 1 GeV 
results in roughly 10% variation in g, 6% variation in 
h whereas the value of g is more volatile and can even 
change sign for large values of fi. 

If we do not include positive parity states' contribu- 
tions in the loops (and naturally fix Ah = 0), we obtain 
a best fit for this coupling g = 0.53. We see that chiral 
loop corrections including positive parity heavy meson 
fields tend to increase the bare g value compared to its 
phenomenological (tree level) value of gi,o. ~ 0.61 [iot . 
while in a theory without these fields, the bare value 
would decrease. 

The fitted value of \h\, is close to its tree level phe- 
nomenological value obtained from the decay widths of 
Dq and D[ mesons (and using the tree level value for 
/ = 130 MeV) hi,o = 0.52. 

Our determined magnitude for g is close to the QCD 
sum rules determination of its dressed value [l^, [4l| , but 
somewhat smaller compared to parity doubling model 
predictions [1, Q. Its sign is also consistent with the 
lattice QCD result of Ref. [S^. Based on this calcula- 
tion we can derive a prediction for the phenomenological 
coupling between the heavy axial and scalar mesons and 
light pseudo-Goldstone bosons Gp*Po^, which is related 
to the bare g coupling as (see e.g. Ref. ji^ ) 



G 



PTPorr 



1 



Pott- 



(13) 



Using our best fitted value for g = —0.06 and excited 
meson masses from TABLE U we predict the absolute 
value of this phenomenological coupling for the case of 
P* = d[°, Po = D*+ and it = ^"t |G^;o^.+,-| = 6.0 
corresponding to an effective tree level coupling value of 

iffr.li = o.i5. 

Next we study the effects of the counterterms on our 
couplings fit. Following the approach of Ref. [1^ we take 
the values of K5, K5, /tig, re^g, 82 + 6^ and 82 + 5'^ en- 
tering our decay modes to be randomly distributed at 
/i ~ 1 GeV in the interval [—1,1]. Near our original fit- 
ted solution, we generate 5000 values of g, \h\ and g by 
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FIG. 3: Effect of the niq and counterterms of the size order 
|k| on the solutions for the couplings g (top), \h\ (middle) and 
g (bottom) as explained in the text. 



minimizing at each counterterm sample. For each so- 
lution also the average absolute value of the randomized 
counterterms (|k|) is computed. We plot the individual 
coupling solution distributions against this counterterm 
size measure in FIG. [31 We see that the inclusion of coun- 
terterms spreads the fitted values of the three couplings. 
From this we can estimate roughly the uncertainty of the 
solutions due to the counterterms to be at the one sigma 
level g = OMt°Z \h\ = 0.47t°;°I and ~g = -0.06t°;°3 
if we assume the counterterms do not exceed values of 
the order 0(1). This result is in a way complementary 
to the study of renormalization scale dependence of our 
couplings' fit. Both are important since although it is 
always possible in principle to trade the counterterms 
contributions for a specific choice of the renormalization 
scale, the latter will be different for different amplitudes 
where the combination of counterterms will be different. 
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CHIRAL EXTRAPOLATION 



integrals of the form 



Next we study the contributions of the additional res- 
onances in the chiral loops to the chiral extrapolations 
employed by lattice QCD studies to run the light meson 
masses from the large values used in the simulations to 
the chiral limit [28, [291 . 



/r(m, A)| A=largc 

,4-D 



[q^ — m^) A 



(1 



q ■ V 



(14) 



As already noted in the previous section, the inclusion 
of heavy excited mesons in the chiral loops introduces 
large scale dependence into the renormalization of the 
coupling constants, due to the large splitting between the 
ground state and excited heavy mesons in the loops. This 
splitting causes the pseudo-Goldstone bosons in the loops 
to carry large momenta. They can be highly virtual or, in 
the cases of PqPtt, P*P*Tr and P^Pqit couplings' renor- 
malizations even real. Such behavior casts doubts on the 
validity of this extended perturbation scheme, as con- 
tributions from higher lying excited heavy meson states 
seem to dominate the loop amplitudes. The problem can 
be explored by analyzing the dimensionally regularized 
loop integral I^'^ {m, A), which can be found in the Ap- 
pendix B. All the other relevant loop integrals can be 
obtained from this one via algebraic manipulation. The 
integral is characterized by two dimensionful scales (m 
and A). In addition chiral perturbation theory requires 
pion momenta (also those integrated over in the loops) 
to be much smaller than the chiral symmetry breaking 
scale A^. The first integral scale m is the mass of the 
pseudo-Goldstone bosons running in the loop. In lattice 
studies, its value is varied and can be taken as large as 
?7i ~ 1 GeV. Within chiral perturbation theory however, 
it is protected by chiral symmetry to be small. On the 
other hand, once A contains the splitting between heavy 
meson states of different parity, it is not protected by 
either heavy quark or chiral symmetries and can be arbi- 
trarily large. Once we attempt to integrate over loop mo- 
menta probing also this scale, we are effectively including 
harder and harder momentum scales in the dimensionally 
regularized expression as this splitting grows. Finally, as 
these approach A^ , the perturbativity and predictability 
of such scheme break down. 

But while the phenomenological couplings' fit seems 
mainly unaffected by such problems (e.g. the results de- 
pend only mildly on the actual value of the mass splitting 
in the range probed), they play a much more profound 
role in the chiral extrapolation. The issue is especially se- 
vere in the case of pions, which due to their small masses 
can even become real inside the loops. This introduces 
uncontrollable final state interactions and causes the chi- 
ral extrapolation to diverge. To mend this situation, we 
attempt on a approximative solution: We expand the in- 
tegrand of If (m, A) over powers of A. We may do this, 
assuming the relevant loop momentum integration region 
lies away from the (v ■ q — A) pole, which is true for chi- 
ral perturbation theory involving soft pseudo-Goldstone 
bosons and for a large enough A. We obtain a sum of 



where the ellipses denote terms of higher order in the 
1/A expansion. This greatly simplified integral has a 
characteristic, that all terms with odd powers of loop 
momenta in the numerator vanish exactly. Thus, the 
first correction to the leading 0(1/A) order truncation 
appears only at 0(1/ A'^). 

The above described procedure is similar to what is 
done in the "method of regions" (see e.g. [1^) when one 
separates out the different momentum scales, appearing 
in problems involving coUincar degrees of freedom. How- 
ever, here we are only interested in the low momentum 
part of the whole integral and assume the high momen- 
tum contributions are properly accounted for in the coun- 
terterms. The leading order term in ([T^ then yields for 
the loop functions 



Ci{x, m) 
C2{x, to) 



TO 

4x 



TO^ , / TO^ 



(15) 



It is important to stress that the relevant ratio for the 
validity of this approach is A/ E.,^ ^ 1 as we are expand- 
ing in powers of loop momentum, not pseudo-Goldstone 
masses. 

This approach can alternatively be understood as the 
expansion around the decoupling limit of the positive 
parity states with the corresponding contributions being 
just a series of local operators with A dependent prefac- 
tors - effective counterterms of a theory with no positive 
parity mesons. Therefore any large deviations of this 
approach from the predictions of a theory without posi- 
tive parity states and with the couplings properly refitted 
would signal the breaking of such expansion and the fact 
that the contributions from "dynamical" positive parity 
states cannot be neglected. 

We compare the above described approximative ap- 
proach with the complete unaltered expressions for the 
coupling renormalization from the previous section. As 
customary we expect the non-analytic chiral log terms to 
dominate the extrapolation, while any analytic depen- 
dence on the pseudo-Goldstone masses can be absorbed 
into the appropriate counterterms. As an example we 
write down the dominating contributions to the chiral 
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log extrapolation of the g coupling 

^PSPbTT' 9 



1 d5f ■ 



m?, dlogm^ 



(4^/)2 



KcKa + KcKb 



1 - 



2 U29 

g - h - 
9 



1 - 



(16) 



In the above expression wc have for the sake of simplicity 
neglected the SU{3) flavor splittings between the heavy 
mesons which are always small compared to A.sh, are 
of higher order in the power counting and vanish in the 
chiral limit. On the other hand one can immediately 
see, that the Ash contributions due to excited heavy 
mesons in the loops seemingly dominate the chiral limit, 
where they diverge. If we instead use the loop integral 



expansion we effectively replace the 1 



terms 



in Eq. (flB]) with m^/4A|^ which on the contrary vanish 
in the chiral limit. Furthermore, these terms then become 
of the form logm^, and thus formally only contribute 
to the next-to-leading chiral log running. Nontheless, 
we keep them explicitly in our further calculations to 
evaluate the size of their contributions. 

In the following quantitative analysis we take the fit- 
ted values of the couplings from the previous section. We 
compare (I) the loop integral expansion, (II) the complete 
leading log extrapolation with positive parity states con- 
tributions, (III) also in the degenerate limit Ash — 0, 
and (O) a chiral extrapolation without the positive par- 
ity states' contributions and with g — 0.53 as determined 
by the coupling fit from the previous section in such a 
scenario. We assume exact SU{2) flavor symmetry and 
parameterize the pseudo-Goldstone masses according to 
the Gell-Mann formulae as [331 



8Ao?ns 



SApms r+1 



(17) 



p 



where 



mu,d 



/nis 



and SXonis/f^ 
0.468 GeV^. Consequently in the chiral extrapolation we 
only vary the scale parameter r - the light quark mass 
with respect to the strange quark mass which is kept fixed 
to its physical value. We again choose fi = I GeV for the 
common renormalization scale. We subtract the coun- 
terterm contributions which scale linearly with r. The 
slope of these can be different between the two scenarios 
and can be inferred from lattice QCD. In Ref. [1^ the g 
and g couplings were calculated on the lattice at different 
r values. However that study used large values of pion 
masses (r ~ 1) where the predispositions for our extrap- 
olation expansion in scenario I are not justified. Also in 
order to use lattice data to extract the couplings from 
such a chiral extrapolation one would in addition need 
lattice results for the h coupling, since it enters in the chi- 
ral logs of the other two couplings considered. Instead we 
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FIG. 4: The g coupling renormalization in D*"*" — > D^tt^. 
Comparison of chiral extrapolation with g = 0.66, \h\ = 
0.47, g = —0.06 and (I) loop integral expansion (black, solid), 
(II) complete log contribution (dark gray, dotted), (III) the 
degenerate limit (light gray, dash-dotted), and (0) without 
positive parity doublets included in the loops [g = 0.53, \h\ = 
0) (red, dashed line) as explained in the text. 



normalize our results for the g coupling renormalization 
in all scenarios at Sr^iiXoms/ f'^ = A'sj^ (corresponding 
in our case to 7'ab = 0.34) to a common albeit arbitrary 
value of Zp, p^^i{ra.h) = 1 and zero slope. We choose this 
scale since it is a branching point between extrapolation 
regions where different approximative limits may apply: 
by region A we denote the low r range, where our 1/A 
expansion (I) is valid for pion loops; conversely we define 
region B above this point, where the full contribution (II) 
and later also the degenerate limit (III) should apply. In 
order to fit our results to lattice data, one would instead 
need to add the (counter) terms constant and linear in 
r to the chiral extrapolation formulae, representing con- 
tributions from s and u, d quark masses. Their values 
could then be inferred together with the values for the 
bare couplings from the combined fit for all the three 
couphngs to the lattice results. 

As an example we again consider the strangeless pro- 
cess —* £)°7r+ in FIG. [H We can see that including 
the complete chiral log contributions from excited states 
in the loops, introduces large (^ 30%) deviations from 
the extrapolation without these states. If one instead 
applies the approximative approach discussed above, the 
deviations diminish considerably. Surprisingly, the de- 
coupling approach produces results which are numeri- 
cally very close to the limiting case of degenerate mul- 
tiplets. One can understand this in terms of the main 
contributions to the extrapolation still eomming from the 



g terms as expected, albeit now with a different fitted 
value for this coupling. In our approximations the h'^ 
terms then actually reduce this discrepancy due to being 
oppositely signed. Such corrections to the running due 
solely to the terms are less then 5%. 

For completeness we also plot the chiral extrapolation 
diagram for the Z'l.^„ , in FIG. (SI The results for 

the chiral extrapolation of the h coupling renormaliza- 
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Calculation scheme g \h\ g 

Leading order 0.61 \40} 0.52 -0.15" 

One-loop without positive parity states 0.53 

One-loop with positive parity states 0.66 0.47 —0.06 

"Effective tree level coupling value derived from one loop calcula- 
tion for the case Dj^ — » DJ^tt^. 

TABLE II: Summary of our results for the effective couplings 
as explained in the text. The listed best-fit values for the one- 
loop calculated bare couplings were obtained by neglecting 
counterterms' contributions at the regularization scale /j. ~ 
1 GeV. 



FIG. 5: Chiral extrapolation of the h coupling renormaliza- 
tion in D^^ — > D'V"'". Comparison of chiral extrapolation 
with g = 0.66, \h\ — 0.47, g = —0.06 and (I) loop integral 
expansion (black, solid), (II) complete log contribution (dark 
gray, dotted), and (III) the degenerate limit (light gray, dash- 
dotted) as explained in the text. 



tion arc very similar to the g coupling case. Here, in all 
approximations, the main contributions to the extrapo- 
lation come from the g^, g^ and mixed gg terms, with 
smaller corrections due to the Ash dependent h'^ terms 
(except for the complete expression scenario III, where 
these terms dominate). 



VI. DISCUSSION AND CONCLUSION 

Within a HHxPT framework, which includes even and 
odd parity heavy meson interactions with light pseu- 
doscalars as pseudoGoldstone bosons, we have calculated 
chiral loop corrections to the strong D*Dtt coupling 
the DqDtt coupling h and the D[DqTt coupling g. 

The calculations have been done by considering first 
the wave functions' rcnormalization. Due to the rather 
large mass splitting between positive in negative par- 
ity states Ash, we find that the perturbation expan- 
sion holds for scales below fi < I GeV, while these 
new strongly scale dependent corrections become large 
at higher rcnormalization scales. From the data for the 
four decay widths (see TABLE [I| we obtain the best- 
fitted values for the bare couplings, which we summarize 
in TABLE |IT1 We are able to determine all the three 
couplings since the contributions proportional to the cou- 
pling g appear indirectly, through the loop corrections. 
One should remember that the quantitatively different 
results of Ref. [2^ appeared before the observation of 
the even parity heavy meson states and in that study 
a combination of strong and radiative decay modes was 
considered in constraining g. 

Since we consider decay modes with the pion in the 
final state, one should not expect sizable contribution 
of the counterterms. Namely, the counterterms which 
appear in our study are proportional to the light quark 
masses, and not to the strange quark mass [1^. Nonethe- 



less the effects of countertcrm contributions in the decay 
modes we analyze, are estimated by making the random 
distribution of the relevant countertcrm couplings. The 
countertcrm contributions of order 0(1) can spread the 
best fitted values of g, \h\ by roughly 15% and g by as 
much as 60%. Similarly, up to 20% shifts in the rcnor- 
malizations scale modify the fitted values for the g and 
\h\ by less than 10% while g may even change sign at high 
rcnormalization scales. Combined with the estimated 
20% uncertainty due to discrepancies in the measured 
excited heavy meson masses, we consider these are the 
dominant sources of error in our determination of the 
couplings. One should keep in mind however that with- 
out better experimental data and/or lattice QCD inputs, 
the phenomenology of strong decays of charmed mesons 
presented above ultimately cannot be considered reliable 
at this stage. 

In the work presented in , the next to leading terms 
(l/mjy) were included in the study of charm meson mass 
spectrum. Due to the very large number of unknown 
couplings the combination of 1/ mn and chiral corrections 
does not seem to be possible for the decay modes we 
consider in this paper. However, the studies of the lattice 
groups [l2, [H, [11] indicate that the 1 / mn corrections do 
not contribute significantly to the their determined values 
of the couplings, and we therefore assume the same to be 
true in our calculations of chiral corrections. 

Due to computational problems associated with the 
chiral limit, lattice QCD studies perform calculations at 
large light quark masses and then employ a chiral ex- 
trapolation m^r ^ of their results. Our analysis of 
such chiral extrapolation of the coupling g shows that 
the full loop contributions of excited charmed mesons 
give sizable effects in modifying the slope and curva- 
ture in the limit m.^. 0. We argue that this is due 
to the inclusion of hard pion momentum scales inside 
chiral loop integrals containing the large mass splitting 
between charmed mesons of opposite parity Ash which 
does not vanish in the chiral limit. If we instead impose 
physically motivated approximations for these contribu- 
tions - we expand them in terms of 1/ Ash - the effects 
reduce mainly to the changes in the determined values 
of the bare couplings, used in the extrapolation, with ex- 
plicit h contributions shrinking to the order of 5%. Con- 
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sequently one can infer on the good convergence of the 
1/AsH expansion. We conclude that chiral loop correc- 
tions in strong charm meson decays can be kept under 
control, give important contributions and are relevant for 
the precise extraction of the strong coupling constants g, 
h and g. 
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APPENDIX A: FEYNMAN RULES 

In deriving the Feynman rules we set the overall heavy 
quark mass scale to a common scale for all processes and 
states under study inducing a mass gap As (^h) terms 
in the propagators of the positive (negative) parity dou- 
blet states due to the relevant residual mass counterterms 
in the Lagrangian Q . The same approach is taken in re- 
gards to the SU (3) flavor symmetry breaking contribu- 
tions, which also induce mass gaps in the heavy meson 
propagators due to relevant 0(rnq) counterterm contri- 
butions in Lagrangian On the other hand, we neglect 
hyper-fine splittings within individual spin-parity heavy 
meson doublets. These arc degenerate at zeroth order 
in the l/mn expansion at which we are working due to 
heavy quark spin symmetry. Following is a list of derived 
Feynman rules used in the calculations in the text. The 
standard -|-zO-prescriptions are implicitly understood in 
the propagators. 



Pa propagator: 



Pa propagator: 
Poa propagator: 
P^a propagator: 
tt' propagator: 

PaPtir' coupling: 
PaPtTr' coupling: 



P0a{v) 



Plaiv) 



■K\k). 

Pa{v) ^P:iv) 



> < 



PoaPihT^'' coupling: 



PLPibT^' coupling: 



-Pa-PofcTT* coupling: 



P„*P*i,7r' coupling: 



■K\k). 

I 

PT(v) fp*^{v) 



■K\k). 

I 

Poa(v) fP^a-") 
^ 

■K'{k). 

ATW few 
— ^ 

Tv'ik). 
Pa{v) f^Pobiv) 



Tv'ik). 

I 



2(fe-u-Ai,-Aa) 

i 

2(fc-u-As-Aa) 

2(fc-D-As-Aa) 



^ab 



— t K VpAab 



f Aab 



2ig livafl h.a \i 
—f-e K VpAab 



2(fe-U-Ajf-Aa) 



APPENDIX B: CHIRAL LOOP INTEGRALS 

Following is a list of loop integral expressions used in 
the text. Our notation follows roug hly that of Ref. [H. 
All expressions already have the infinite part of the inte- 
grals subtracted as explained in the text. 



d^g 



1 



/r(m,A)=M 



(4-13) 



(27r)^ (g2 - m2) 



IGtt- 



2 1 tilt. 

;m log — 



(27r)^ (q2 _ m2) (w • g - A) 167r2 



Ci ( -,m) q''" + C2 (-,m | V^v' 
ml \m 



d^q 



{2t:)D {q2 - m^){v q - A) 16n^ 



C — ,m 



A 



/3^^(m,Ai,A2) = /i 



{4-D) 



d^g 



Ai - As 



{2Tr)D (q2 -m^){vq- Ai){v ■ q ~ A^) 
[/r(m,Ai)-/r(m,A2)], 



(Bl) 
(B2) 
(B3) 



(B4) 
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where 



/3^''(m,A,A) = — /rKA). 
In the text we then make use of the following expressions 



C(x, to) 



n3 



9 



(B5) 



(B6) 



with 



Ci{x,m) = — 



-12.T + 10.T^ + {9x - 6a:^) log ( — ) - 12x{x - 1)F ( - 



C2(x,m) = C(x,to) — Ci(x,to), 



C'l 2 {x,y,m) = —[Ci^2{y,m) - Ci_2(a;, to)], 

m X — y 



(B7) 
(B8) 

(B9) 



Ci oix, to) — C\ 2 (a;, X, to) ~ C'l 2(2;, to). 

TO dx 



The function F{x) was calculated in Ref. [2£ 







5 (.T + V 


-1). 






f - tan^i 


(v'l-:rO] 







\x\ < 1. 
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